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^ ■ Abstract 

O. 

Organic heterojunction solar cehs are analyzed within a minimal model that includes the essential 



physical features of such systems. The dynamical properties of this model, calculated using a master 
equation approach, account for the qualitative behavior of such systems. The model yields explicit 
results for current-voltage behavior as well as performance characteristics expressed in terms of the 
thermodynamic efficiency as well as the power conversion efficiency at maximum power, making it 
possible to evaluate the optimal setup for this device model. 
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I. INTRODUCTION 

The limited supply of today's main energy sources, such as oil and coal, will force us to 
rely increasingly on renewable energy sources. Besides wind and water power, energy con- 
version based on photovoltaic (PV) devices has received much attention as such a source. 
Promising systems for next generation devices are organic photovoltaic solar cells (OPV) 
because of their potential for low-cost processing. Of particular interest are polymer-based 
hetero junctions consisting of a blend of electron-donor (D) and electron-acceptor (A) mate- 
rial (for reviews, see Refs. [l|, |2| , and references therein). A prominent acceptor material is 
the buckminsterfullerene (Cgo) 3|-[6|. The quest to improve the energy conversion efficiency 
of such systems is the focus of intensive current research. 

To evaluate and subsequently improve the efficiency of organic photovoltaic cells it is cru- 
cial to understand the underlying energy conversion processes and how material properties 
affect their overall performance. Widely accepted is the multistep generation process that 
starts with photon absorption by the donor (often a polymer) yielding an exciton (bounded 
electron- hole pair). The generated exciton diffuses to the D-A interface, where it dissociates 
into free charge carriers, which are later transported to the electrodes. The D-A interface 
should be constructed to favor energetically fast and efficient electron transfer leading to 
exciton dissociation. 

The dynamics of electron transfer at the D-A interface is of crucial importance for the 
performance of hetero junction solar cells as measured by their efficiency. In considering 
this issue, one may address the thermodynamic efficiency 77* 7|, |8| by considering the solar 
cell as a heat engine operating between a hot and a cold reservoir with temperatures Tg 
( "sun temperature" , representing the incident radiation) and T (temperature of the chemical 
environment), respectively. Alternatively, the conversion efficiency 77 [9] associated with the 
maximum power point in the current- voltage {J -U) characteristic is of interest as a realistic 
performance measure. Establishing the relationship between system properties that affect 
exciton dissociation at the D-A interface and the cell efficiency is a major goal of the ongoing 
research. Within this effort, it is useful to consider simple model systems for which one can 
obtain explicit relationships between system structure and characteristic parameters and its 
performance measures. In this article, we describe and analyze such a model system. 

Our model consists of coupled donor and acceptor molecules, each described as a two level 
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FIG. 1: A minimal model of an organic heteroj unction PV cell. The system consists of one donor 
(e.g. a suitable polymer) and one acceptor, e.g. fullerene sites, each characterized by their HOMO 
and LUMO levels. 

(highest occupied molecular level, HOMO, and lowest unoccupied molecular level, LUMO) 
system, situated between two electrodes. As such, it is an extension of a simpler model 



^ 



but with an important additional 



recently analyzed in a similar context by Rutten et al. 
feature - the existence of an heterojunction characterized by energetic parameters - Coulomb 
interaction and donor-acceptor LUMO-LUMO gap, that were identified as important driving 
factors in the operation of such systems. The system dynamics associated with this model is 
described by a kinetics scheme derived using a lattice gas approach. 
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12l | similar in spirit 



to previous work ISMlTJ] that uses a master equation approach to analyze cell dynamics. 

We show that an effective mechanism for both exciton pair formation and dissociation can 
be captured by a minimal model of this type, which can by used as a framework for discussing 
the current- and power-voltage curves and the cell efficiency. In particular, the model leads 
directly to the predictions of an optimal interface energy gap Ae (usually associated with the 
energy difference between the lowest unoccupied molecular levels (LUMO's) of the D and A 
molecules) for these efficiency measures. This can be compared with the performance of an 
ideal device, in which nonradiative losses are absent. In this ideal case, the thermodynamic 
efficiency is found to decrease monotonously, whereas the power conversion efficiency still 
goes through a maximum, with increasing Ae. 



4 5 

FIG. 2: The six accessible microstates of the model system considered. The probability to find the 
system in state j (j = 0, ..., 5) is denoted Pj. 

II. MODEL AND COMPUTATIONAL DETAILS 

A. Model 

The PV cell model considered here is a one-dimensional lattice gas with two "effective" 
sites / = D, A. One site represents the donor (D, e. g., polymer-based material) and the 
second acts as acceptor (A, e. g., fullerene-based material). Each of the sites is represented as 
a two-state system with energy levels (£:di)^d2) and (£aij^a2) corresponding to the (HOMO, 
LUMO) levels of the donor and acceptor species, respectively. In what follows we use the 
notation AEi = ei2 — en {I = D, A), for the energy differences that represent the donor 
and acceptor band gaps, and refer to Ae = e^a — ^a2 as the interface or donor-acceptor 
LUMO-LUMO gap. The microstates of the system can be specified by the set of occupation 
numbers n = (nm, n^a, fij^^, UAa), where n^ = or 1 {I = D,A; j = 1, 2) if the corresponding 
level is vacant or occupied by an electron. 

To assign further realistic contents to this model, we introduce the restrictions nDi^D2 = 
and tiai = 1. The first of these restrictions implies that the donor can be in the ground or 
excited state, or, following charge separation, in positively ionized state, but excludes its 
doubly occupied (negatively charged) state. The excited donor state (noi, Uoa) = (0, 1) rep- 
resents the exciton formed as result of light absorption. The second condition implies that the 



TABLE I: System states and their occupations with corresponding energies. 



STATE 


OCCUPATION 

(rioi, flD2) ^Al) ^A2) 


ENERGY 





(0,0,1,0) 


£0 = £ai 


1 


(1,0,1,0) 


£1 = £di + £ai 


2 


(0,1,1,0) 


£2 = £d2 + £ai 


3 


(0,0,1,1) 


£3 = £ai + £a2 + Vc + V^ 

= Eai + £a2 + V'c 


4 


(1,0,1,1) 


£4 = £di + £ai + £a2 + Vc 
= £di + £ai + £a2 


5 


(0,1,1,1) 


£5 = £d2 + £ai + £a2 + Vc 
= £d2 + £ai + £a2 



acceptor can be in either its ground (n^i, ^^2) = (1, 0) or its negative ion (^^1,^1^2) = (1, 1) 
states. Therefore, the system is characterized by six states with respect to the occupations 
{riui, nj^2, n^i, ^^2), that we denote by the integers 0, ..., 5 as shown in Table [H 

In the expressions for the states energies, Vc > is the Coulombic repulsion between two 
electrons on the acceptor, whereas V^ > is the Coulombic energy cost to move an electron 
away from the hole remaining on the donor. In general, we expect that Vc > V^. The 
sum Vc + V^ can be thought of as the exciton binding energy in this model: It is the total 
Coulomb energy cost for dissociating the exciton on the donor by moving an electron to the 
acceptor. It is convenient to redefine ej^2 = ^a2 + Vc so that the energies Ej {j = 0, ..., 5) are 
determined by the five energy parameters £01, £02, £^ai, £'a2, and V^. 

At the left (donor, say) and right (acceptor) ends of the system, the device is connected 
with two electrodes represented by free electron reservoirs at chemical potentials (Ik {K = 
L, R). This is a highly simplified picture that disregards electron transport within the donor 
and acceptor phases. We have opted to make this simplification to focus on the important 
step of interfacial exciton dissociation, but future more realistic treatments should take these 
components of the overall dynamics into consideration. As sketched in Fig. [H we use the 
common picture by which the left reservoir is assumed to exchange electrons only with the 
HOMO level of the donor, whereas the right lead exchange electrons with the upper level 



of the acceptor. Direct electronic interaction between system and reservoirs is not explicitly 
taken into account, but it is implicit both in the states of molecular species adjacent to 
metal electrodes and in the kinetic charge transfer rates. 

Next, we construct the kinetic scheme for the time evolution of the average occupations 
Pj (j = 0, .., 5) associated with this level structure (see Fig. |2]). In writing these equations, 
we make the simplifying assumption that electron exchange between molecules and metals 
involve only metal electrons at the electrochemical potentials /Xl and /Xr of the left and right 
leads (corresponding to a bias potential U = (/Xr — fii^)/\e\ where e is the electron charge. 
We also disregard possible environmental relaxation dynamics due to polarization effects 
associated with the formation of transient charged molecular species. Generalizing this 
dynamical picture to take such processes into account (e.g. by considering time dependent 
site energies [l8[) is another important subject for future work. The kinetics process in our 
scheme then corresponds to the following processes: 

(a) Electron transfer between levels Dl and the left electrode and between A2 and the right 
electrode, with rates determined by the corresponding molecular energies, molecules- 
leads coupling, electrochemical potentials in the leads, and the environmental temper- 
ature, T. 

(b) Electron transfer between donor and acceptor, determined by the coupling between 
them, the corresponding state energies, and the temperature, T. 

(c) Light induced electron excitation (rate kg) and relaxation (rate kg) between donor levels 
Dl and D2. These rates are modeled as thermal rates determined by the corresponding 
state energies and radiative coupling, and the "sun temperature" T^. 

(d) Radiationless (thermal) electron transitions between level Dl and D2 with rates k^,. 
(excitation) and fcnr (relaxation), determined as in (c) except that relevant coupling is 
vibronic in origin and the temperature involved is the environmental temperature T. 



Explicitly, the transition rates kj/j = kj/^j from state j to j' are given by 

^01 = ^34 = K = ^'l/I^l) , (1) 

kw = hs = k^ = z/l/(xl) , (2) 

kl2 = "^45 = "^S + "-nr 

= I/s(l + ^s(a;s)) + t'nrll + nnr(a;nr)) , (3) 

/i;21 = hi = ks + /^nr = l^sns{Xs) + t'nr^nrla^nr) , (4) 

^32 = ^DA = I'uaK^va) 5 (5) 

^23 = W = l^UAfiXuA) , (6) 

A;o3 = ku = k25 = K = ^bRxj,) , (7) 

kso = ^41 = k52 = k^ = z/r/(xr) . (8) 

In these equations /(x) = 1/ [exp(x) + 1] and /(x) = 1 — /(x), Xl = (^di ~ /^l)/^b7', 
a^R = {^ A2 — f^n) / ksT , and Xda = (£^d2 — £'a2 — V^O/.'^bT'. The rate coefficients z/ are determined 
by the corresponding couplings. Following Ref. ?]], the transition rates between the HOMO 
and LUMO in the donor phase are assigned by using the boson population factors ng = 
1/ [exp(xs) — 1] and n^^ = 1/ [exp(xnr) — 1] with scaled energies Xg = AEj^/kj^T^ and Xnr = 
AEjy/k-gT. /cl = koi = A;34 is the rate to move an electron from level Dl to the left electrode, 
/cr = A;o3 = ^14 = ^25 is similarly the rate from A2 to the right electrode, /cl = kiQ = ^43 and 
k^ = ^30 = ^41 = ^52 are the corresponding reverse rates. The transition rate from D2 to A2 
is given by k^A = ^32, whereas the opposite rate is /cda = ^23- Also, we denote by ks and fcnr 
the radiative and thermal excitation rates in the donor and by kg and fcnr the corresponding 
relaxation rates, so ks + A^nr = ^21 = ^54 and fcg + k^^ = ki2 = ^45. 

B. Computational Details 

The master equation with the rates Eqs. (l)-(8) that describe the time evolution of the 
probabilities Pj = P^n^i, rioa, "^ai, '^a2) (j = 0, ..., 5) to be in the six possible states thus 



reads 



dPojt) 

dt 
dPiJt) 

dt 



dP^jt) 
dt 



dPsjt) 
dt 



dPi{t) 
dt 



koiPiit) + ko3P3it 

hoPo{t) + kuPiit 

- {koi + k^i + /cai 

k2lPl{t) + k23P3{t 

- {ki2 + h2 + h2 

ksoPoit) + k32P2{t 

- (^03 + ^23 + kij, 

k^sP^it) + k^,Pl{t 
-{ku + ki4 + k^i 



- {kio + kso) Poit) 
+ kuP2{t) 

+ k25P5{t) 
P2{t) , 

+ k^^P^{t) 
Psit) , 

+ A;45A(t) 
P^it) , 



and normalization implies that 



P5(t) = l-$^P,(t). 
i=0 

In terms of these probabilities electron currents can be expressed as 

J^{t) = k,{Po + Ps)-k{Pl+P4) 

Ut) = K{Ps + P4 + P5) - kn{Po + Pi + P2) 

Jsit) = ksiPl + P^) - ksiP2 + P5) 
Jnrit) = k^,{P, + P4) - ~knr{P2 + P5) 

JuAyt) = kTj^P2 — kiy^Ps . 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 
(16) 
(17) 
(18) 
(19) 



Jh {Jr) is the current entering (leaving) the molecular system from (to) the electrodes, Js 
and Jnr are, respectively, the light induced and nonradiative transition currents between the 
HOMO and the LUMO in the donor phase, and Joa is the average current between the 
donor and acceptor species. Below we will focus on the steady state magnitude, Jl = Jr. = 
Jda = Js + Jnr = J) of these currents. 

The time evolution and steady state currents associated with this kinetic model can 
be evaluated exactly, however such exact solution becomes costly for larger, more realistic 
system that takes into account also transport within the donor and acceptor phases. We 



therefore advance also an approximate mean field treatment and show that it can provide 
a good approximation to the exact analysis and can be used for realistic multilevel systems 
at a relatively low computational cost. 

To this end, we introduce the averaged site occupations Poi{t) = {n^j^) = Pi{t) + Pi{t), 
p^,{t) = {n^.,) = P2{t) + P5(t), and p^^{t) = {uj,,) = Ps{t) + P^{t) + P^{t) (see also the 
illustration in Fig. |2]). Therefore, the currents (15)-(19) can be rewritten as follows 

Ji.{t) = K{pm - Pm) - hPoi , (20) 

Js{t) = ksPo, - ksPu2 , (22) 

where p = 1 — p. The treatment of the average current JoA{t) between the donor and 
acceptor species is more difficult. In terms of occupation number JoA{t) it is 

■JoAit) = koA{nu2 (1 - ^A2)) - /^Da((1 " ^Dl " ^02) ^A2) • (24) 

Neglecting fluctuations, i. e., {nj^^rij^^) ~ (^di)(^a2) = PdiPa2 and {nj^zUj^^) ~ (^d2)(^a2) = 
Pd2Pa2, leads to the mean fleld expression of the current 

■^DA V") ~ "^DaPd2Pa2 ~ "^Da(Pdi ^ Pd2)Pa2 ■ (25 j 

As a consequence, we arrive at the following mean fleld rate equations for the averaged site 
occupations 

dpoi 



dt 
dpo2 

dt 
dpA2 



Mt) - Js{t) - Jnr(t) , (26) 

Mt) + Jnrit) - JZ{t) , (27) 

JZ{t)-Ut). (28) 



dt 

The full kinetics in this approximation is obtained by solving Eqs. (26)- (28) together with 
Eqs. (20)- (23), and (25) self consistently. Note that the number of coupled equations solved 
in this scheme grows linearly with the number of N of single electron states, while in the 
exact approach this number is essentially the number of molecular states ~ 2^. 

To illustrate the nature of the kinetics that results from these rate processes, we choose 
a set of physically reasonable parameters: In the calculation discussed below the following 



< 

c 



^ 







8 




— J..,.t 






I/,„ = O.OlyDA \ \ 


6.4 




— P^^. 






\ ■; 


4.8 




-^MF 

^MF 








3.2 








1.6 




(a) ^^- 


^ 




V 





^ 











16 

- 12.8 

9.6 

6.4 

3.2 





c 



ttn 



0.2 0.4 0.6 0.1 

t/(V) 



1 1.2 1.4 






7.U 




8 

6.4 
4.8 
3.2 
1.6 




^ 


"exact n ni V 
P ^ i'n. = O.OlfDA V 

-fexact J Y 
<^ exact 1 ^ ^iT^, V 

(b) ^'^ 1 



16 

12.8 

^ 9.6 

6.4 

3.2 







0.2 0.4 0.6 0.1 



1 1.2 1.4 



FIG. 3: Current (J, left vertical axis) and power (P = ?7J, right vertical axis) in the PV cell, 
plotted against the voltage bias U (see text for parameters). Plot (a) compares results from the 
exact calculation to those obtained from mean field approximation. Plot (b) compares the same 
exact results (i^nr = O.OIz^da) to those obtained for the ideal cell for which z/nr = 0. We assume 
that z^DA = lO^^s"^ 



choice was used: /11 = 0.0 eV, /in = /11 + |e|f/, ^di = - 
l^=_0.25eV, and Vc' = 0.15 eV. Thus, AEj, = £d2-£di 



0.1 eV, e^. 



1.4eV, e, 



0.9 eV, 



1.5 eV [19( and 1^0 + l^c = 0.4 eV 

20|, |2l|. For the temperatures we take T = 300K and Tg = 6000K. The kinetic rates 

are set to z^l = z/r = i/g = ^'nr = O.OIz/qa and z/qa = 10^^s~\ describing a system with 

efficient and fast donor-to-acceptor electron transfer (which can occur on the ps timescale 



22]), and moderate radiationless losses, as would be used in such applications. Finally, 



note that in the particular example employed here we have considered a situation where the 
imposed potential bias falls between the acceptor species and the right electrode. It should 
be emphasized that although the results shown in Figs. 3 and 4 are based on these choices, 
the qualitative behavior discussed below holds for a wide range of these parameters. 
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III. RESULTS AND DISCUSSIONS 

Figure [3] shows results for the stationary donor — )■ acceptor current obtained from both the 
exact solution from Eqs. (9)- (19) and the mean field approximation, Eqs. (20)- (28). Figure 
12^1' compares the results of the mean field approach to their exact counterpart, showing that 
the former provides a good approximation to the exact behavior. Figure [3]d compares the 
performance of a junction characterized by the prescribed parameters to the corresponding 
ideal junction in which nonradiative losses are absent (z/nr = 0). In both cases, the current 
is constant until the electrochemical potential on the right electrode comes within ~ fegT of 
the acceptor level A2, and decreases sharply after it exceeds this level. Consequently, the 
generated power, P{U) = U J{U) goes through a sharp maximum, Pmax = ^max <^(^max) in 
that voltage region. 

Next, consider the efficiency. The maximal power conversion efficiency is defined by 

■""^^ P. ■ (2'" 

where Pa, is the incident radiant power, a constant independent of the process undergone by 
the system. The thermodynamic efficiency at maximum power is given by 

if -'max L'max ^l^f-^maxj ^max / -. ^nrl^^-^maxj \ /'Q^^ 

where Qs is the net energy absorbed per unit time from the radiation field, that is, heat 
absorbed from the hot reservoir, Qg = AEj^Jg. We note in passing that all processes 
undergone by our system are accompanied by well defined energy changes: The heat fiuxes 
associated with the electron exchange processes between system and leads, are, in the present 
model Ql = (^di — f^hjJh and Qr = (/Ur — eA2)<^R, Qm = '^E-^Jj^j. is the net heat generation 
per unit time related to nonradiative relaxation processes, and Qda = — Ae: J^a is the heat 
fiux associated with the electron transfer at the D-A interface. Energy conservation implies 
that the overall cell power is the sum of these fiuxes, i. e., P = Ql + Q^ + Qs + Qm + Qda = 
if'ii- l^i.)J = UJ. 

Of central importance is the dependence of the efficiency on the interface LUMO-LUMO 
energy gap Ae. One may intuitively expect to find that an optimal value of this parameter 
exists: A finite Ae is needed to drive the charge separation process but a larger Ae implies 
that more energy may be lost to unproductive processes. Equations (29) and (30) quantify 
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this phenomenon, which is illustrated in Fig. |H Figure H^ shows the power conversion 
efficiency calculated for the chosen parameter set using Eq. fl29|) . For comparison, the 
corresponding result for the ideal cell {u^j- = 0) is also shown. Both are seen to go through 
a maximum as functions of the interfacial LUMO-LUMO gap Ae. The thermodynamic 
efficiency, Eq. f l30p . displayed against Ae in Fig. |Ud, shows a similar pronounced maximum, 
however the ideal thermodynamic efficiency is a monotonously decreasing function of Ae. 
Note that the ideal cell efficiency is an upper bound to the actual efficiency, however, for 
this finite power operation, it is below the Carnot efficiency, rjc = 1 — T/T^ = 0.95. 

Finally, consider again the performance of the mean field approximation relative to the 
exact solution. As seen from Figs. 3 and 4, the mean field treatment provides a good approx- 
imation that closely follows the behavior of exact solutions. The influence of correlations 
is particularly seen in Fig. HJd, where we find small but noticeable differences between the 
mean field and the exact curves for Ae beyond the maximum at ~ 0.4 eV. 

IV. CONCLUSIONS 

In summary, the effort to increase and optimizing the efficiency of heteroj unction OPVs 
necessarily involve many structural and energetics system parameters. The above consider- 
ations focus on what is arguably the most important generic issue - the optimization of the 
interfacial LUMO-LUMO gap that compensate between the need to overcome the exciton 
binding energy and the required minimization of losses. For the present simplified model 
and our choice of parameters we find the most efficient setup for interface gap energies 
somewhat above 0.4 eV. More important is the fact that the present model with future gen- 
eralizations that should include transport in the donor and acceptor phases and polaronic 
relaxation following redistribution of charge densities provides a framework for analyzing 
such efficiency measures. The mean field approach introduced here provides a reliable ap- 
proximation that can be used for fast evaluation of more complex model systems and will be 
useful in extending these studies to realistic OPV cell models. A more ambitious task would 
be to generalize the concept advanced in the article to the quantum regime. Consideration 
of quantum effects coherence (see, for example, Ref. 
of fundamental limits to photovoltaic efficiency. 
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23l ) may prove useful in the discussion 
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FIG. 4: (a) The power conversion efficiency t] and (b) thermodynamic efficiency rf evaluated 
at maximum power and displayed as functions of the interface gap Ae. Mean field results are 
compared with the exact solutions of the underlying master equation. Parameters are the same as 
in Fig. [3] with the difference that eA2 runs from 1.4 eV to — 0.4eV and /i^ is determined by -Pmax- 
For the power, P = P^ incident on the cell we assume 1.0 nW j24i] . 
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